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Abstract 



Three-dimensional scalar electrodynamics, with a local U(l) gauge symmetry, is believed to 
^ ' be dual to a scalar theory with a global U(l) symmetry, near the phase transition point. The 

' conjectured duality leads to definite predictions for the scaling exponents of the gauge theory 

transition in the type II region, and allows thus to be scrutinized empirically. We review 
these predictions, and carry out numerical lattice Monte Carlo measurements to test them: a 
number of exponents, characterising the two phases as well as the transition point, are found 
to agree with expectations, supporting the conjecture. We explain why some others, like the 
exponent characterising the photon correlation length, appear to disagree with expectations, 
unless very large system sizes and the extreme vicinity of the transition point are considered. 
Finally, we remark that in the type I region the duality implies an interesting quantitative 
relationship between a magnetic flux tube and a 2-dimensional non-topological soliton. 
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1. Introduction 



Dualities constitute one of the very few analytic tools available for studying non-perturbative 
properties of systems with many degrees of freedom. They have found applications in widely 
different physical settings, ranging from spin models to quantum field and string theories. 
A duality transformation typically maps topological defects to fundamental fields, and vice 
versa, and can therefore translate a non-perturbative problem to a solvable perturbative one, 
either in the same or in a different theory. 

In some spin models, the duality transformation can be carried out explicitly, mapping the 
partition function of one theory to that of another (see, e.g., Refs. The duality is 

therefore a mathematical identity, and one can see exactly how the parameters and observ- 
ables of the two theories relate to each other. Interacting continuum field theories, on the 
other hand, are generally fairly difficult to treat exactly. One reason is that they contain 
fluctuations on many different length scales, which do not decouple from each other, and this 
leads, among other things, to ultraviolet divergences. The duality may then be approximate 
rather than exact, and precisely valid only in a certain limit in the parameter space. 

One way to limit the effect of ultraviolet divergences is to lower the dimensionality, and for 
instance the sine-Gordon and Thirring models in 1-1-1 dimensions have famously been shown 
to be dual to each other [SUHI- In higher dimensions, a way to control ultraviolet divergences 
is to consider supersymmetric theories, where they are weak or even absent. For instance, the 
M = 2 SU(A''c) super- Yang- Mills theory exhibits a Montonen-Olive duality |H1E1 between 
electric charges and magnetic monopoles. Recently there has also been a great deal of interest 
in dualities between the large- A'^c limit of four-dimensional = 4 super- Yang-Mills theory 
and string theory in five-dimensional anti-de Sitter background |10j . 

Clearly, it would be interesting to extend a quantitative understanding of dualities to- 
wards dimensions closer to the physical 3-1-1, or to less supersymmetric theories, in order 
to approach, for instance, a consolidation of the 't Hooft-Mandelstam dual superconductor 
picture of confinement in QCD (for a recent review, see Ref. ^J). So far, however, rather 
few examples are available in these directions. In this paper, we study one of them, a duality 
between two non- supersymmetric continuum field theories, namely scalar electrodynamics 
(SQED) and complex scalar field theory (SFT), in three Euclidean dimensions [H], jl2j-j2Uj. 
The duality maps the Coulomb and Higgs phases of SQED to the broken and symmetric 
phases of SFT, respectively. The microscopic details of the theories are not dual to each 
other, but as one approaches the transition point, the duality should become a better and 
better approximation at long distances. In principle, the duality becomes exact at the tran- 
sition point (in the so-called type II region), and also describes correctly the approach to the 
transition point. 

The nature of the duality becomes more transparent if one considers the theories in 2-1-1 
dimensional Minkowski space. Then the duality maps the fundamental fields of one theory 
to vortices of the other: the vortex lines can be understood as world lines of particles in the 
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dual theory. One can easily make some elementary observations that hint towards the duality. 
First, both the Coulomb phase of SQED and the broken symmetry phase of SFT have one 
massless degree of freedom, namely the photon and the Goldstone mode, respectively |18j . 
while in the other phase these degrees of freedom go over into two degenerate massive modes, 
in both cases. Furthermore, both the interaction between vortices in SFT and the Coulomb 
interaction between electric charges in SQED have the same logarithmically confining form. 
In contrast, the Abrikosov-Nielsen-Olesen vortices in the type II region of SQED and the fun- 
damental particles of SFT with a positive quartic coupling, have an exponentially decreasing 
repulsive Yukawa interaction. 

A full duality would, of course, be a stronger statement than simply a counting of the 
degrees of freedom, and mean that the two theories have identical dynamics as quantum 
systems. If all the operators allowed by the symmetries were kept in the two theories and 
one knew exactly the mapping between their parameters, the duality would predict that for 
any observable there is a corresponding observable in the dual theory with the same value. 
For instance, the mass spectra of the two theories should be identical. 

In practice, however, one of the theories in question is truncated by dropping an infinite 
series of high-dimensional operators, and the exact mapping between the parameters related 
to the operators kept is not known, because there is no supersymmetry to cancel radiative 
corrections. Therefore, the best way to see the duality is to calculate universal quantities such 
as critical exponents, which describe the scaling properties of the system as the transition 
point is approached. Many critical exponents of SFT are known to a high accuracy |21j . 
because the theory is in the same universality class as the three-dimensional XY model. The 
duality predicts that each of these exponents has a dual counterpart in SQED, which should 
have exactly the same value. ^ 

In this paper, we invoke previously developed numerical techniques to study topological 
defects j22[ I2,'S| and certain two-point functions j24| l2f)j in SQED, to measure a number of crit- 
ical exponents with lattice Monte Carlo simulations in the type II region, and compare them 
with predictions following from the duality conjecture. The purpose is to demonstrate that 
these techniques yield results precise enough to serve as very non-trivial checks of this con- 
jecture. The paper is organised as follows. In Sec. El we formulate the duality transformation 
and show how it relates the critical exponents of the two theories. We describe our numer- 
ical simulations and present their results in Sec. El In Sec. 01 we elaborate briefly on some 
interesting qualitative manifestations of the duality conjecture in the case of a macroscopic 
external magnetic field, both in the type I and in the type II regions of SQED, although we 
have not carried out any new simulations for this situation. Finally, we discuss our findings 
and present our conclusions in Sec. El 

^To be precise, the SQED observables we consider are sensitive only to the SFT exponent u-^y as well as 
a certain anomalous dimension rj, while for instance the exponents /3xYi7XY related to the response of SFT 
to an external magnetic field, do not play a role in the following 0. 
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2. Duality 



2.1. Basic setup 

The gauge theory under investigation, three-dimensional (3d) SQED, is formally defined by 

CsQEB = \Fii + iDk(t>nDk(t>)+m^cP*<p + \{^*(t>f , (2.1) 
^SQED = J VAkV4>V<t>* exp(^- I^Csqed) , (2.2) 

where F^i = dkAi — diAk, Dk = dk + ieAk, k,l = 1, 2, 3, repeated indices are assumed to be 
summed over, and = Jd'^x. This theory is super-renormalisable, with the only divergences 
appearing in the parameter m?. Unless otherwise stated, we assume to denote the bare 
parameter, and the theory to be regulated, for the moment, with dimensional regularisation. 
Note that all the arguments that follow are non-perturbative in nature, and there is thus no 
need for gauge fixing in Eq. H2.2|) . 

While the duality relation can be expressed explicitly in a certain deformation of SQED, 
namely when the theory is discretised and the so-called London limit (Aa — > cxd, where a 
is the lattice spacing) is subsequently taken [3], mi-^Z], its status is less clear in the full 
continuum theory of Eqs. (|2.1() . H2.2|) . where A is finite and o — > 0. In fact, the phase diagram 
of SQED is even qualitatively different from that of the London limit, where the transition 
is of the second order. In SQED the transition is of the first order if the ratio A/e^ is 
small |26| (the type I region), and is believed to be of the second order for a large A/e^ above 
the Bogomolny point (the type II region). This belief is based, apart from the apparent 
vicinity of the London limit, on 3d renormalisation group studies [231-1101) previous lattice 
simulations ^S], as well as other arguments ^^-^ 

To now formulate the duality conjecture, let us choose a specific set of observables to act 
as probes for the properties of the system. Defining a magnetic field through 

= ^eijkFjk{x) , (2.3) 

we will mostly consider objects related to the two-point function (see also Refs. [241 \'25\ ) 

Ckiix - y) = (Bk{x)Bi{y)) . (2.4) 

To compute this kind of objects, let us generalise Eq. (|2.2|) and define 

^SQED[^i] = j VAkVcpVf exp[- l^i^CsQEB - Hi{x)B,{x))] , (2.5) 

where Hi{x) is an arbitrary source function. The object in Eq. (|2.5() is in principle well defined 
for a complex Hi{x), but in the following we restrict ourselves to a stripe in the complex plane 

■^As is well known, the renormalisation group in 4 — e dimensions, with e <C 1, continues to predict a 
first-order transition even at large A/e^ [2^1 E3- 
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where Hi(x) consists of a real constant part (or zero-mode) , representing a physical external 
magnetic field strength, and a purely imaginary space-dependent part, used as a probe to 
define various Green's functions. (If one wishes, the zero-mode could also be considered to 
be purely imaginary to start with, and analytically continued to real values only in the end, 
at least as long as the volume is finite.) The reason for this choice is that it turns out to 
simplify the form of the corresponding SFT, without posing any restrictions on the physical 
observables that we can address. The correlator of Eq. (|2.4|) is now obtained by taking the 
second functional derivative of In ^^sqed [-f^i] ) and setting Hi = afterwards. 

The basic observation behind the duality is that SQED allows to define an exactly conserved 
"charge", the magnetic flux through some surface. This conservation law can be expressed 
in a local form by noting that, identically, 

d^Biix) = . (2.6) 

Thus, there should exist some global symmetry for which Bi is the Noether current . 
This global symmetry should be broken in the "normal" Coulomb phase where the U(l) 
gauge symmetry is restored, the massless photon representing the Goldstone boson ^Hj. In 
the "superconducting" Higgs phase where the U(l) gauge symmetry is broken, on the other 
hand, the global symmetry should be restored, and the particle number of SFT, representing 
the number of Abrikosov-Nielsen-Olesen vortices in SQED, locally conserved. 

If we were to stick to the Coulomb phase only, which has a single exactly massless degree 
of freedom, then the fact that an effective description exists can be made rather rigorous, 
following the general arguments for effective theories Once we approach the phase 

transition point, however, other (non-Goldstone) degrees of freedom also become massless, 
and the situation is less clear. 

The assertion of the duality conjecture is that ZsQED[Hi] equals the partition function of 
a scalar field theory (SFT), with explicit global U(l) symmetry and a mass parameter 
which is (almost) the reverse of m^, 

= Co — cim^ + ... . (2.7) 

Here ci is dimensionless and positive, and cq is constrained by the requirement that (with a 
given regularisation) , = cq — cim^ + ..., where m^, are the values of the mass parameters 
at the respective transition points. The quartic coupling A of SFT is also constrained: 

~X = do + di\ + ... , (2.8) 

where di is dimensionless and positive, and do is constrained by the requirement that Ac = 
do + diXc + where Ac is the "tricritical" value separating the type I and type II regions 
in SQED JH]) while Ac = is the value at which scalar particles turn from attractive to 
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repulsive in SFT.^ The conjecture then reads that 



^sqed[^^] =^sft[^.] , (2.9) 

where the form in which Hi appears in the partition function Z^ft remains to be determined. 

In order to make progress, let us note that the symmetry underlying the duality can be 
considered to be a local one, if we assign a transformation law also to the source field Hi (cf., 
e.g., the analogous procedure in the case of chiral symmetry in QCD; Ref. [HH] and references 
therein). Indeed, we can allow, in general, the source term Hi to change by a local total 
derivative, since this leaves 2sqed invariant, due to Eq. (|2.6|) . provided that the boundary 
integral emerging vanishes. Let us inspect this issue in some more detail. 

To understand the significance of the boundary term, let us place the system in a finite 
box of volume V, and choose boundary conditions such that all physical fields, like Bi{x), 
are periodic, in order to preserve translational invariance. In accordance with our choice 
for the analytic structure of Hi(x), the transformation property of Hi has now to be purely 
imaginary, 

Hi{x) H[{x) = Hi{x) - idioix) , (2.10) 

where a{x) is an arbitrary real function, defined, up to a so far unfixed proportionality 
constant, to be the generator of the local U(l) symmetry transformation. The partition 
function defined in Eq. I|2.5)) has in fact a larger symmetry than U(l), corresponding to a 
complex function a(x), but as we shall see, it is only transformations of the type in Eq. (|2.1()j) 
which match the properties of the U(l) symmetry on the side of SET. Now, since a{x) is not 
directly a physical field, we can in general consider the possibility that its boundary conditions 
are not strictly periodic, but are non-periodic in some direction, by an amount which we call 
A for the moment. This leads, in what one might call "large gauge transformations", to 



^SQED[i^-] = ZsQEDlHi] exp A y"d2s . B 



(2.11) 



where the integral is over that boundary of the box at which a{x) is discontinuous by A. We 
observe that if we choose A = me, where m is an integer, and the usual flux quantisation 
condition e / d^s.B = 27rn, with n an integer, then the system is indeed fully invariant within 
each "topological" sector, characterised by n. (The whole partition function contains a sum 
over all the sectors). In the following we set m to its lowest non-trivial value, m = 1, so that 
the discontinuity is A = e. 

On the SET side, then, we assume the symmetry transformation generated by a{x) to 
operate explicitly on the dual field variable, (j): 

4)^4)' = e'^'^4> 4>* 4>*' = e-'^'^4>* , (2.12) 



^If SFT is used for describing weakly interacting atomic Bose-Einstein condensates, then the scalar self- 
coupling is usually written as A = where a is the s-wave scattering length and m the atom mass; 
the transition from attractive to repulsive corresponds to the s-wave scattering length changing sign. 
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where we have introduced a proportionahty constant e. We now see, first of all, that ea is 
defined only modulo 27r. Therefore, boundary conditions can also only be imposed modulo 
27r for ea, and we obtain a relation to the constant A introduced above, 

^=1 = ^. (2.13) 

A e 

Furthermore, since ZsYT:[Hi\ must be invariant under the local transformation of Eqs. (|2.10j) . 
H2.12P like ZsQED[-ffi] is, we can finally fix its structure: 

CsF^{H,) = \ ZFl + [{dk - ~eHkW][{dk + eHk)$\ + rh^*^ + ~X{4>*4>f + ... , (2.14) 

2sFT[i^^] = y'pWexp[-^£sFT(^0] , (2.15) 

where Fki = dkHi - diHk- 

Inspecting Eq. ()2.14|) . we immediately obtain a physical interpretation for the duality. 
Choosing a real constant and coordinates so that ^ only in one direction, say X3, 
which we may rename to be "imaginary time", we note that Eqs. (|2.14|) . (|2.15j) represent 
just the Euclidean path integral expression (cf., e.g., Ref. [lOj) for a complex SET with the 
chemical potential fi = eH^ related to the conserved current 

jk = 2Im[4>*dk4)] = 4>idk4>2 - <p2dk4>i , (2.16) 

where we have written (p in terms of real components ascj) = ((^1+2(^2) /V^- The corresponding 
total particle number / d^s thus represents in SET the same object as to which jj, = eH^ 
couples in SQED, namely the integer-valued conserved flux (e/27r) / d^s • B, where d^s is the 
spatial volume element. 

As we have written SET in Eq. 1)2. 14() . it only contains a few terms. In principle, how- 
ever, Eq. 1)2. 14() should include an infinite series of (gauge invariant) higher order operators, 
starting with ~ {(j)*(f))'^ . In the type II region (A > 0), they are just at most marginal, and 
do not modify any of the critical exponents at the transition point. In other words, their 
contributions are suppressed by some power of M/A, where M denotes the dynamical mass 
scales inside the truncated action of Eq. 1)2.14(1 , while A ~ is a "confinement" scale related 
to excitations within 3d SQED that remain non-critical at the transition point. In the type 
I region (A < 0), on the other hand, the higher order operators are important. In fact, A < 
together with a positive (</)*(^)^-term provides just the usual prototype for a first order phase 
transition, as is the case in the type I region. In most of the discussion that follows, we 
consider the type II region, and can thus ignore the higher order operators. These arguments 
also implicitly assume that the external magnetic field eHk is "small", i.e. at most of the 
same order of magnitude as the dynamical mass scale M. Note that for a constant eH^, 
terms like -F^^ vanish. 

At the end of the day, the external source field Hi is often set to zero, and in that case, 
the Lagrangian vCsft of Eq. p.l4|) simplifies to the standard one, 

CsFT = dk4>*dk4> + rh^*4> + ~X{4>*^f + - • (2.17) 
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The general form in Eq. (|2.14|) is still important, however, because typical predictions of the 
duality conjecture follow by taking functional derivatives of the equality in Eq. ()2.9|) . and 
setting Hi ^ afterwards. In particular, the second functional derivative produces for our 
main probe, Eq. (|2.4j) . 



Cki{x - y) 



(5^ lnZsQED[Hi 



6Hk{x)6Hi{y) J 

= z{d 6ki - dkdi)6{x -y) + 2eHki5{x - yW4>{x)) - e^{jk{x)ji{y)) + ... , (2.18) 

where the expectation values on the right-hand side are evaluated with the Lagrangian of 
Eq. (|2.17|) . Both sides of the relation in Eq. (|2.18|) are transverse. The constant Z is seen 
to contribute to "contact" terms only, ~ 6{x — y), but it is significant if the corresponding 
susceptibility (integral over all space of the two-point correlation function) is considered, 
which thus is not determined by the dual theory alone. Non-contact terms, on the other 
hand, are fully predicted in terms of the parameters of Eq. 1)2. 17(1 . 

To conclude, we should reiterate that in certain limits, for instance when SQED is re- 
placed with a "frozen superconductor" , or integer valued gauge theory, relations of the type 
in Eq. (|2.18|) can be made exact, for any values of the parameters (the only one being the 
inverse temperature (3 in that case) jHlESl- In our case, on the other hand, the relations be- 
tween the parameters are largely open, and the parameters even get renormalised differently 
in the two theories. It is only the infrared properties of the theory, that is correlations of 
the type in Eq. (|2.18j) at non-zero distances and close to the transition point, which can be 
related to each other, by suitably tuning the coefficients Cj in Eq. 1)2. 7() and di in Eq. (|2.8() . 

2.2. General structure of the photon two-point correlator 

Defining now the Fourier transform 

Bk{p) ^ I e'P-^B,{x) , (2.19) 



the object we will mostly consider is the Fourier transform of Eq. (|2.4() . 

CUp) = y{Bki-p)Biip)) = l^e^P<y~^\Bk{x)Bi{y)) ^ (^^ - ^)g(p) , (2.20) 

where V is the volume. We choose to measure the correlator in the xs-direction and use a 
momentum transverse to that direction, p = p (for instance, p = 27rnl/L, where 1 is the unit 
vector in the xi-direction, L is the linear extent of the system, and n is an integer), so that 

C33(p) = G(p) . (2.21) 

The general structure of G(p) is 
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where, on the tree-level in SQED, the self-energy is just a constant, S(p) = my, with mv 
the inverse of the vector (or photon) correlation length. More generally, we expect that close 
to the critical point, 

S(p)^m| + A|p|2-'' + 0(|p|^), (2.23) 

where A is some (dimensionful) constant, r] is the anomalous dimension, and 5 > 2 — r]. 

Taking a Fourier transform of Eq. (|2.18|) . e*^'^^"^-* (...), we obtain now a prediction based 
on duality for the observable of our interest, G{p): 

G(p) = -Zp^ + 2e\r^)-e^y{j3i-p)j3ip)) + ..-. (2.24) 

The non-trivial object here, (|j3(— p)j3(p)^, has properties determined by SFT alone, and 
leads thus to a definite structure of G(p), including corrections to scaling, etc. Formulating 
Eq. 1)2. 24() in such an explicit form is, as far as we know, a new result in the present context. 

In principle, Eq. 1)2. 24(1 could now be used to obtain a direct prediction for the object G(p), 
by inserting the numerically measured properties of SFT on the right-hand side. In practice, 
this is not quite straightforward, because the relations of the parameters are not fixed by the 
conjecture. What can be done, however, is to use Eq. 1)2.24(1 to obtain predictions for various 
critical exponents, since such predictions are parameter-free. SFT is known to be in the same 
universality class as the three-dimensional XY model, many critical exponents of which are 
known numerically very well [2^. Therefore, we will use them as a benchmark. 

2.3. Critical exponents 

While SQED in the type II region allows to define a large number of critical exponents, we will 
in this paper restrict only to a few of them, related in one way or the other to the "magnetic" 
properties of the theory, such as in Eq. (|2.4|) . The reason is that the corresponding observ- 
ables can be measured with controllable statistical and systematic errors, thanks to various 
numerical techniques introduced in Refs. j22j-|25j. and that analytic predictions for these 
observables, following from the duality conjecture (through Eq. (|2.24l) '). are unambiguous. 
We start by considering the symmetric phase. 

Symmetric phase: magnetic permeability. Magnetic permeability can be defined by 
considering a constant source = H, and the response of the magnetic field strength 
B = V-^J^B3{x) to H, 

dB 1 9^ 

^^dH = VdiP ^"^^Q™ = i^o^(P) • ^2-2^) 

In the free Abelian theory we would have H = B and x = !> as can be seen by setting E — > 
in Eq. (|2.22() . or directly by starting from the free energy F{B) = | B^V. When interactions 



8 



are taken into account, 1-loop perturbation theory at large gives 



e2 



+ ... , (2.26) 



247ry^^?(^ 

where m?{p,) is the renormalised mass parameter in dimensional regularisation (in the MS 
scheme, for concreteness) . Thus, x decreases as we approach the transition point m?{j2) ~ 0. 

To understand the critical behaviour of Xi note that in the broken phase of SFT, the 
correlator {jk{—p)jiip)) is dominated by the massless Goldstone mode, so that {jk{—p)3i{p)) oc 
{(t)*(j)) PkPi/p^, and consequently 

Iim(j3(-P)i3(p)) =0. (2.27) 

p->0 

Thus the duality in Eq. (|2.24|) predicts that 

X = 2e2(r0). (2.28) 

Essentially the same argument was presented recently by Son [20]: he considered an effec- 
tive theory for the Goldstone mode alone, (/>(x) = (l/\/2)/ exp(ia;(a;)), and thus obtained 
X = e^/^, where is the helicity modulus, or stiffness, or decay parameter, related to the 
Goldstone mode, and equals 2{(f)*(p) in our notation. Note that we have implicitly assumed 
the use of continuum (dimensional) regularisation in writing down the duality relation, and 
Eq. (j2.28j) should also be understood to be correspondingly regularised. 

There are various ways to derive the critical scaling exponent of For instance, 

given that the action contains S'sft ~ / d^x dk(l)*dk(j>, we may dimensionally expect jJQ that 
{(f)*(j)) ~ l/\x\ ~ 1/^sFT ~ 1''"^-^^) where we have introduced the distance from the critical 
point, r, by 

2 2 

m — mr. 

-T^ = r . (2.29) 



Other arguments leading to the same result can be put forward (Ref. [211] and references 
therein), and will also be met below. Thus, if x{'t) ~ kT'^ foi' 't ~^ O^i we obtain [20] 

i^x = ^XY , (2.30) 

where i^xY is the exponent of the correlation length ^sft in SET. 

To test this critical behaviour we can devise, following |21], a finite-size scaling procedure for 
measuring the exponent. Close to the critical point, for box sizes larger than the correlation 
length ^ of non-Goldstone modes, the system behaves as if it were almost in infinite volume, 
so that (pmin = 27rl/L) 

S(Pmin) = X"'(r)pL„ + O(Pmin) • (2-31) 

Eor box sizes smaller than ^ (but still large compared with the microscopic scales), on the 
other hand, the system behaves as if it were already at the critical point, 

S(p^i„) = ^p^T^^ + 0{pIJ . (2.32) 
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Therefore, 

PLi'^"^(Pmin) ~ P;nTnS(pmin) ~ I ^ ' L<^ ' ^^'^^^ 

The function has to be continuous at L ~ ^ ~ |r|~'^XY; thus, using Eq. ()2.3U|) . r] = \. The 
functional form close to the critical point should be universal, and according to Eq. ()2.33p 
for ?y = 1, only dependent on \t\L^/'^^: 

PminG-l(p^in) ~ f{\T\L^'''^) , (2.34) 

with f(x) = A, X <^ 1, and f{x) ~ x~'^^, x ^ 1. Subsequently at the critical point, |r| — > 0, 
it has a fixed point value, A. Moreover, if we take a derivative with respect to r at the critical 
point, we obtain 

~ L^/"^ . (2.35) 

r=0 

Eq. (|2.35|) will be used below to measure i'^, as was done in Ref. for a related model. 



u 
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Transition point: anomalous dimension. In the previous paragraph, we already found a 
specific value for the anomalous dimension rj, rj = 1, assuming that the magnetic permeability 
scales with the same exponent as the correlation length, as argued in Eq. (|2.3U|) . Let us 
now show that the same value for the anomalous dimension can be obtained without any 
assumption for the behaviour of x{t)- The reasoning could then be reversed, to provide yet 
more evidence for the scaling of x{t) according to Eq. (|2.3U|) . 

The anomalous dimension related to G(p) can be found in a particularly simple way with 
an argument similar to one by Son [20]. Let us consider ZsFT[Hi] at the transition point. 
Factoring out the term multiplied by Z, the system should be "conformally invariant" at 
the critical point, or have no scales. For dimensional reasons, the quadratic part of In^^gpT; 
which is a function of eHi only, must thus have the structure 

ln^SFT[^*] oc / (6ki - ^) [eHk{-p)][eHi{p)]\p\ , (2.36) 

where = Jd?p/ {2t[Y . Using this to compute G(p) through the Fourier transform of the 
second functional derivative of InZgpT, hke in Eqs. 1)2. 18() . ()2.2U|) . 1)2. 24(1 . and comparing 
with Eq. fTT^ for m| = 0, it follows that r/ = 1. 

The structure of Eq. 1)2. 36(1 can also be found with a perturbative 1-loop computation. On 
the side of SQED, the computation was carried out in Ref. 02] in the Coulomb phase, with 
the result ^ 

S(p) = ^ e^lpl , for p2 » m\jx) . (2.37) 
lb 

Higher loop corrections are non-vanishing, however, and even diverge at the critical point 
m?{p,) ~ 0, whereby there is really no hard prediction for the critical exponent. On the 

^Wilson-type renormalisation group studies in SQED lead to the same functional behaviour |28l I29| . 
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other hand, one could also carry out the computation on the side of SFT, as we just did, 
and possibly combine with the Wilson renormalisation group there: this line of reasoning can 
in principle also be used to understand that ry = 1 (see, e.g., Ref. [19 J. For a lattice study 
within a certain version of the dual theory, again leading to = 1 , see Ref. fi3\ . 

Broken symmetry phase: inverse vector propagator at zero momentum. The 

quantity we will measure in the broken symmetry phase is the infrared limit of the inverse of 
the vector propagator, 

lim p^G"^(p) = lim S(p) = m| . (2.38) 

p— >0 p^O 

Our aim is to determine the exponent associated with m|, m| ~ \Tp^. To relate 7s to more 
physical quantities, note that according to Eq. 1)2. 23() . for small p and the structure of 
the inverse propagator should be 

p2G-i(p) -ml + AIpp-" . (2.39) 

Therefore the inverse of the physical vector correlation length, or the "photon mass" my, 
which is defined by the position of the singularity in G(p), or zero in p^G~^(p), scales as 

mv = IPpolel ~ (m|)^ ~ ^ Uv = . (2.40) 

2 — r] 

Thus, a determination of 7s combined with the known rj amounts to a determination of uy- 
Conversely, given a prediction for z/y, as follows in the next paragraph, we have a prediction 
for 7s ^1 , which will be tested below. 

Broken symmetry phase: vector correlation length. As just mentioned, the vec- 
tor correlation length is determined by the position of the singularity in G(p). According 
to Eq. (|2.24)) . it is determined in SFT by the singularity structure in the current-current cor- 
relator in the symmetric phase. Since the current represents a two-particle state (~ (pid(j)2, cf. 
Eq. (|2.16l) '). it is natural to expect that the singularity is placed at the two-particle threshold, 
my = 2M, where M denotes the inverse of the scalar correlation length; this is certainly the 
behaviour obtained with a 1-loop computation in SFT. Therefore, we expect that 

mv = 2M ~ 2|r|''XY , (2.41) 

i.e., vy = i^XY |291 119j . Returning now back to Eq. 1)2.40(1 and inserting ij = 1, we also find 
that 121 

7s = uxY . (2.42) 
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Figure 1: Left: Examples of reweighted susceptibilities x( I </'P)- Right: Determination of the 
infinite volume critical point yc (triangle) from the positions of the susceptibility maxima. 

Broken symmetry phase: vortex tension. Consider finally the vortex tension. Follow- 
ing again an elegant argument by Son j^Hj, let us consider the effect of a constant = H on 
the gauge theory side. Because of the Meissner effect, the system does not respond (^sqed 
does not change) until H > Hd = eT/27r, where T is the tension of an infinitely long vortex. 
On the SFT side, on the other hand, a constant eH corresponds to a relativistic chemical 
potential /i in a (2+l)-dimensional theory, fi = eH. Therefore, the system does not respond 
until fi > M, where M is the particle mass. We thus obtain M = eHd = T, i.e., that the 
vortex tension should again scale with the same exponent i/xY as M does: if T ~ \t\'^'^ , then 

i^T = iyxY. (2.43) 



3. Simulations 

3.1. Discretised action 

We discretise the action in Eq. (|2.1|) in a standard way, by replacing 

1 



1 



Fh(x) — ^ afc(x) + a/(x + a/c) - afc(x + a/) - a/(x) = — ^aM(x) , (3.1) 

Dk4> - exp(iafe(x))(/)(x + aA;) - 0(x) , (3.2) 
a L J 

/d^x ^ ^a3, (3.3) 
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Figure 2: Left: The behaviour of G(pmin) at a few representative volumes. The filled 
circles and the solid curve correspond to the infinite volume limit. Right: The function 
Pmm^^^(Pmin)- The solid curve corresponds to the infinite volume limit. 

where a is the lattice spacing, k,l are unit vectors, and ai(x) = aeAi{'x.).^ The bare mass 
parameter is written in a form which guarantees that the MS renormalised mass parameter 
m?{fi), evaluated at the scale fl = e'^, remains finite in the continuum limit 

2 _ 2/ 2^ . 2 ^^^^ 3.175911535625 
m — m [e ) — (e + 2A) — 



27ra 

1 



167r2 



Ae^ + 8Ae^ - SA^) ( In + 0.08849 ) - 1.1068e^ + 4.6358Ae2 



(3.4) 



The couplings e^,A, on the other hand, do not require renormalisation in three dimensions 
(for 0(a)-corrections, see Ref. |17|). Once we also rewrite (p = ecj), the action becomes 
dimensionless, parameterised only by 

The continuum limit corresponds to Pg ^ oo. Instead of Eq. 1)2. 29(1 . we can now write 

T = y-yc- (3.6) 

The fields that are updated are the real variables afe(x) and the complex ones </>(x). For 
details concerning the update algorithm employed in this work, we refer to Refs. j22ll23j . Once 
the system is put on the lattice, one also has to impose boundary conditions. In the following 
paragraphs we employ periodic boundary conditions for all the fields. This implies that the 
net winding (or flux) through the lattice in any configuration, w = e /d^ s-B, vanishes, rather 

^Note that in Eq. I|3.1^ we use the non-compact formulation for the U(l) gauge field. We do this to avoid 
topological artifacts, monopoles, which appear in the so-called compact formulation |44) . unless e^a <^ 1, and 
would make it difficult to reach large volumes in physical units, i.e. V ^ |45) . 
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Figure 3: Left: Examples of infinite volume extrapolations for Pjjjfjj^(Pmin) in the symmetric 
phase. We fit to the 1-loop perturbative expression, with the mass of the scalar particle as 
a free parameter. The extrapolation range increases rapidly while approaching the critical 
point. Right: The corresponding behavior for G{0) = limp^^.^^o[l + PmL^(Pmin)]^^- The 
dashed line represents the perturbative prediction in Eq. 1)2. 26() . with m?{fL) — > (y — Uc)^'^- 



than fluctuates as it in principle should in the canonical ensemble of Eq. 1)2. 5|) . At the critical 
point this choice affects some quantities, like amplitude ratios, but it is not expected to affect 
the critical exponents on which we concentrate here. 

The observables are discretised in a straightforward way. In particular, denoting again 
p = 27rnl/L, the quantity G(p) in Eq. ()2.21)) is measured as 

G(p) ^ /?G E e*P"(ai2(0)ai2(x)) . (3.7) 

X 

In the figures we always plot the momentum in the form in which it appears in discrete space, 

2 . a|p| 2 . a|pmin| 2 . vr 

PA = - sm — - , pA.min = - sm = - sm — , (3.8) 

a 1 a 2 a Jy 

where p was assumed to point along the xi-axis, and the lattice size was denoted by L = Na. 
We often use furthermore lattice units, a = 1. 

3.2. Simulation parameters 

All the simulations in this paper have been carried out with a single lattice spacing, f3G = 1. 
The reason is that since we are interested in universal critical behaviour and are using the 
non-compact formulation, there is no need for a continuum extrapolation, as long as there 
is no phase transition in between the /3g used and /3g = oo; this indeed is the case. The 
scalar coupling is chosen to lie comfortably in the type II region (i.e., beyond the Bogomolny 
point, X = 1/2), X = 16. We have also performed some simulations at x = 2, which still lies 
in the type II region confirming the qualitative pattern observed at x = 16 but with 
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Figure 4: The derivative —dy[pmmG ^(Pmin)] at the critical point, as a function of the 
volume, and a determination of the exponent of magnetic permeability via Eq. ()2.35|) . 

significantly less resolution. Volumes have been chosen in the range 6^. ..52^, and for each 
parameter value we collect from ~ 10^ to ~ 2 x 10^ sweeps. Different parameter values are 
joined together with Ferrenberg-Swendsen multihistogram reweig hting [IH]. 

3.3. Location of the critical point 

The first task is to locate the critical point, ye- While there are many possibilities for doing 
this, all are in principle equivalent in the limit y — > oo, and we choose here to use the location 
of the maximum of the "specific heat", or the susceptibility related to 

xM')^N'[{{Wf)-{Wf] , (3.9) 

where |(/>p = 4'*4'- Our data at a few representative volumes are shown in Fig.^left), 

and an infinite volume extrapolation based on the positions of the susceptibility maxima in 
Fig. fright). The extrapolation has been carried out with the finite-size scaling ansatz 

y,{N = oo) = y,{N) + + C2j^ + ... , (3.10) 

where the exponents have been fixed to their SFT values, u = uxy = 0.67155(27), uj = loxy = 
0.79(2) HU. We find 

yc = -17.749(5) . (3.11) 

This value will be frequently referred to below, in the form yc ~ —17.75. 

The critical point is of course clearly visible also in the observable we are actually interested 
in, G(p). In Fig. 121 we show the structure of G(pnim) (left) and Pmin^~"'^(Pmm) (right), 
with Pmin = 2tt1/L. Both become order parameters in the infinite volume limit Pmin — ^ 0, 
vanishing on one side of yc- 
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Figure 5: Left: The photon self-energy S(p) at the critical point. Circles denote momenta 
|Pmin| = 27r/L at different volumes; triangles denote larger momenta. The solid curve is a 
linear fit determining the anomalous dimension rj. Right: An alternative determination of r] 
as well as of the fixed-point coefficient g* (cf. Fig. inj, based on Eq. (|2.33|) . 



3.4. Symmetric phase: magnetic permeability 

The magnetic permeability, as defined in Eq. 1)2. 25() , could in principle be determined directly 
from an infinite volume (pmin ^ 0) extrapolation of the data in Fig. [2tleft). Such an ex- 
trapolation cannot be carried out in practice, however, because the function approaches its 
infinite volume limit very slowly close to the critical point. This is illustrated in Fig. |21^1eft) 
at two selected y- values well inside the Coulomb phase, y > Uc- The corresponding extrapo- 
lated values are shown in Fig. ISfright) (without any estimate of the systematic uncertainties 
introduced by the extrapolation), but are reliably extracted only so far above yc that no 
meaningful fit for a critical exponent can be carried out. 

Fortunately, we can determine the exponent directly, by carrying out a finite-size scaling 
study exactly at y = yc, employing Eq. 1)2. 35(1 . In order to measure the derivative, we 
construct explicitly the corresponding operator, 

-ry°" = c^i^G = -^[(i;e*-ai,(0)a,,(x)i:#-(y)#(y)) 

5;]e*"a,2(O)a,2(x))(^0-(y),^(y))] .(3.12) 
X y 

The data are shown in Fig. ^left), and a fit to the filled circles produces a value 

= 0.659(26) . (3.13) 
This is well consistent with the SFT value vxY = 0.67, as predicted by Eq. (|2.30jl . 
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Figure 6: Left: Determination of the fixed point value of the photon self-energy from the 
function S(pmin)/|pmin|5 by a comparison of two volumes. Right: An extrapolation to infinite 
volume and a comparison with the determination of g* in Fig. [S] (open triangle). 

3.5. Transition point: anomalous dimension 

Our next task is to determine the anomalous dimension. In Fig. Efleft), we do this directly 
through Eq. with i;(p) extracted from the measured G(p) (Eq. ^J^) via Eq. 

A fit to the data at Pmin produces 

rj = 0.98(6) , (3.14) 

in perfect agreement with the prediction following from Eq. H2.36() . 

An alternative determination of rj can be obtained by using Eq. ()2.33|) . In Fig. [SJ^right), 
we show S(pinin)p^7n V = Vci observing that a fixed-point value in the infinite- volume 
(Pmin — ^ 0) limit is only obtained with ry ~ 1. We have also estimated the fixed-point value 
g* of S(pinin)/Pmin, Corresponding to A/e^ in the notation of Eq. (|2.33j) .^ 

The fact that the function S(pjnin)/Pmin approaches a fixed-point value, allows also for an 
alternative determination of g* . Indeed, following the standard procedure, we can estimate 
the infinite-volume limit by comparing data for S(pjnin)/Pmin at two lattice sizes, N and 2A'^, 
for various N. The procedure is illustrated in Fig. IH^left), and produces values for g* (cf. 
Fig. Efright)) consistent with but more precise than in Fig. [SJ^right): in the limit A'^ — > oo, a 
linear fit yields g* w 0.13(1). This value will find applications in Sec. 13.71 

3.6. Broken symmetry phase: inverse vector propagator at zero momentum 

We next consider the parameter m^, defined in Eq. (|2.38|) . The infinite- volume extrapola- 
tion of p^^;^G~^(pinin) is illustrated in Fig. [Tfleft), and a fit as a function of y is shown in 

®The fixed-point value g* , unlike critical exponents, is possibly sensitive to the boundary conditions used. 
Let us reiterate that our simulations at this point correspond to a vanishing winding (or flux), w = 0. 
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Figure 7: Left: Infinite volume extrapolations for p^G^^(p) at two values of y in the broken 
symmetry phase. Circles denote momenta |pmin| = 2tt/L at different volumes; triangles 
denote larger momenta. Right: The scaling behavior for the infinite volume extrapolated 
p^G^^(p). The inlay shows the dependence of the critical exponent on the value of Uc- 

Fig.Etright). We find 

7s = 0.651(33) , (3.15) 

in perfect agreement with Eq. (|2.42)) . 

3.7. Broken symmetry phase: vector correlation length 

Let us then discuss the vector correlation length. As reported earlier on ^Hl OSji direct 
measurements of the vector correlation function struggle to show critical behaviour according 
to the exponent vxY as predicted by Eq. 1)2. 41(1 . producing rather the exponent i/xy/2. In the 
light of the data presented in this paper, however, such a behaviour is well understandable. 
Indeed, close to the critical point the form of the vector propagator is as shown in Eqs. 1)2.22(1 . 
()2.23j) . with T] = 1. We have just determined that A ~ g*e'^ ~ 0.13e^. This means that in 
order for true critical behaviour to show up in infinite volume, we would need to require 

ml < 0.136^7711/ , or 777V < 0.136^ . (3.16) 

This is the case, however, only extremely close to the transition point (cf. Fig. 2 in Ref. |35j). 
Moreover, even when 777y ~ 0, finite-volume corrections to the asymptotic form of the photon 
correlator are important, unless 

p^iin <C 0.13e^pmm , (3.17) 



which at Pg = 1 transforms to 



N > « 48 . (3.18) 

0.13 ^ ^ 
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Figure 8: The photon correlation function at the critical point, with transverse momentum 
Pmini compared with a Fourier transform of + ^l^*!] (solid; cf. Eqs. H2.22() . ()2.23|) ). as 
well as with Fourier transforms of (dashed), l/|p| (dotted), in each case with an overall 
constant chosen so that data points are matched at r/aN = 0.5. 



It is not easy to satisfy this inequality in practice, however. Away from the transition point 
my ~ 0, or in volumes smaller than Eq. 1)3. 18() . on the other hand, the free term will 
dominate the inverse vector propagator in Eq. (|2.22j) . and solving for the pole position with 
the free part alone, we recover the exponent 

.f-*^™)^l^^ = i.XY. (3.19) 

For completeness, the highly non-trivial form of the vector correlation function is illustrated 
in Fig. El at y = yc- It is seen clearly how the free term p^jj-, and the linear term in $](piuin) 
are both needed in order to reproduce the data points, even at volumes as large as 52^. 

If the mapping between the parameters of SQED and SET were known, one would not 
need to rely on the asymptotic critical behaviour to demonstrate the duality between the two 
theories, but one could compare the data directly, staying away from the transition point and 
keeping the volume finite. Without the mapping such a comparison is not available on a quan- 
titative level, but we may still note that the data obtained for the vector correlation length in 
the "frozen superconductor" model, where the duality is exact, show a behaviour very similar 
to what we have argued for here and observed in Refs. jHSlESl, i-e. ^'xy/2 < z^^ffective ^ u^y^ at 
least apart from the extreme vicinity of the transition point j2H] • This provides further qual- 
itative support for the duality conjecture, through similar corrections to asymptotic scaling 
in both theories. 
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T/e^ T/e^ v=0 ,671 (38) 




10 20 30 N "'^ "'^^ y-17.4 



Figure 9: Left: An example of an infinite volume extrapolation (triangle) for the tension 
T/e^ in the broken symmetry phase. The ansatz is as explained in Ref. j^, and is only 
applicable in large volumes (solid curve). Right: The scaling behavior for the infinite volume 
extrapolation. The inlay shows the dependence of the critical exponent on the value of ye- 

3.8. Broken symmetry phase: vortex tension 

We finally consider the vortex tension, T7 Unlike for the observables so far, its determination 
requires that we depart from strictly periodic boundary conditions for all the fields. In fact, 
the Monte Carlo determination of T proceeds via a computation of the free energy difference 
IuZsqedI'U^ = 0) — lnZsQED(w^ = 27r), where the winding number w = e /d^s ■ 'B = 27r 
corresponds to the presence of a single vortex through the lattice. Parametrizing w = 2ttz, 
we can write lnZsQED('U^ = 27r) = In^sQED(O) + ^Qdz{W{z)), where {W{z)) represents 
a specific expectation value, which can be measured by Monte Carlo methods. For more 
details, we refer the reader to our earlier work |22| I23j . 

In the earlier work mentioned, the integral ^q(1z(W {z)) was approximated by a sum over 
independent and consecutive Monte Carlo measurements at a discrete set Zj, i = l,...,n^, 
where the parameter had typically values = 11, ...,21. For the current work, we have 
implemented "parallel tempering" in the parameters Zj, along the lines of Ref. [IHI) in an 
attempt to sample intermediate configurations with non-integer winding numbers < z < 
1 more efficiently. Using again = 21, we find a sizable error reduction, as compared 
with our earlier approach. This is illustrated by the single data point with large error bars 
in Fig. I^right), at y ~ —17.9, obtained with the unimproved algorithm. We have not 
investigated the precise mechanism behind this algorithmic improvement in great detail, 
however. 

In Fig. in^left) we show an example of the infinite volume extrapolation of the tension in 

''We use the same notation for the vortex tension and the temperature, but there should be no danger of 
confusion, since the two are never discussed in the same context. 
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the broken phase, and in Fig. I^fright) a fit to the extrapolated values. The fit yields 



i^T = 0.671(38) , (3.20) 

in perfect agreement with Eq. (|2.43|) . 

4. Macroscopic external magnetic fields 

Up to this point, we have mostly considered the case that the magnetic field is set to zero. 
Hi = 0, after taking derivatives of the partition function (cf. Eq. (|2.18|) ): the main exception 
was the discussion of the vortex tension, leading to Eq. (|2.43|) and continued in Sec. 13.81 In 
this section, we wish to make some further qualitative remarks on the case that the external 
magnetic field is not set to zero but kept finite. 

To begin with, let us reiterate the duality relation for this situation, in rather explicit 
form. We choose again coordinates such that the constant magnetic field is pointing in the 
third direction, = H, and consider the system to live in a finite box (or hypertorus), 
now of extent L in the third direction. Then Zsqed[H] represents the classical Euclidean 
partition function for 3- dimensional SQED, with the external parameters H,L. At the same 
time, the corresponding ^sft[^]> Eqs. (|2.14|) . (|2.15j) . is just the imaginary time quantum 
partition function for a 2-dimensional SET, in the presence of a finite chemical potential 
fi = eH (written in a relativistic form) and a finite temperature T (fixed by the extent of 
the imaginary time direction). Thus, the relations between H, L on one side and fi, T on the 
other, read 

e he , , 

Taking the box cubic and sending L ^ oo, corresponds to computing the quantum partition 
function for the 2-dimensional SET at zero temperature. While the relations in Eq. ()4.1|) 
have been expressed in a canonical ensemble with respect to H, /i, it should be clear that the 
correspondence remains true also in a microcanonical ensemble (i.e., fixed magnetic flux in 
SQED / fixed particle number in SET), with a proper choice of boundary conditions. 

Implications for the type I region. In the type I region of SQED, vortices attract each 
other. Therefore, if we are in the broken symmetry phase and force a magnetic flux through 
the system, the vortices form a flux tube that penetrates through the Meissner phase. What 
is the analogue for this in SET? The type I region corresponds to A < (and the presence 
of further stabilising terms), and as is well known, in this case the system admits non- 
topological soliton solutions [EOIHI) or droplets of Bose liquid, to be intuitively thought of 
as non-dispersive bound states of particles held together by the attractive interaction. We 
thus observe a perfect analogy, with the role of flux quanta in SQED played by the attractive 
particles in SET. Note that even though SET appears in a (2-1-1 dimensional imaginary time) 



21 



relativistic form in Eqs. ()2.14|) . H2.15() . the same solutions appear there even if we approach 
the non-relativistic hmit |511 I52j . a situation also familiar from actual experiments with 
soliton-like structures in superfluid Helium (see, e.g., Ref. [SSI) atomic Bose-Einstein 
condensates (see, e.g., Ref. '5l|). 

Implications for the type II region. In the type II region of SQED, on the other hand, 
vortices repel each other. Therefore, if we place the system in a magnetic field, they form (at 
least on the mean field level) an Abrikosov vortex lattice. It has been a long-standing issue to 
study whether, within SQED, the vortex lattice indeed exists in a strict sense for low enough 
magnetic fields, and then melts to a vortex liquid possibly through a first order transition, as 
observed experimentally [SS], or always appears in a liquid state due to fiuctuations, so that 
the transition observed in experiments would be a manifestation of some physics beyond the 
pure SQED (see, e.g., Ref. [Sni)- We now see that on the side of SET, this corresponds to 
whether a 2-dimensional dilute system of atoms with repulsive interactions, forms a lattice, 
or melts due to quantum and thermal fluctuations. In principle, this issue could be studied 
with SET, more easily than directly with SQED as attempted in Ref. since the system 
has fewer dynamical length scales. In practice, though, the inclusion of a chemical potential 
in SET makes the action complex for a generic configuration (this is just the generic "sign 
problem" for // 7^ 0), rendering importance sampling ineffective. 

5. Conclusions 

In this paper, we have used numerical lattice Monte Carlo simulations to measure a number of 
critical exponents in three-dimensional scalar electrodynamics (SQED) in the type II region, 
or A/e^ > 1. We have also reiterated which exponents of the dual scalar field theory (SET) 
they should correspond to. Our measurements agree well with the SET values, known to high 
accuracy from previous studies of the XY model. This provides strong "empirical" evidence 
for the existence of the conjectured duality between the two theories. 

We have also elaborated on the qualitative implications of the duality for the case of 
macroscopic external magnetic fields, both in the type I and in the type II regions. Even 
though we have no new simulations to report in this regime, we find it remarkable that for all 
the known phenomena observed in SQED somewhat below the phase transition temperature, 
one indeed finds a perfect counterpart in 2-dimensional quantum SET, but at a different 
temperature, as determined by the physical geometry of the SQED sample. 

The duality predicts that the only important degrees of freedom near the transition point 
in SQED are vortex lines, and the quantities we have measured are sensitive to their proper- 
ties. The agreement with the dual theory therefore also demonstrates that the methods we 
have used in these simulations give an accurate description of the dynamics of topological 
defects. Similar methods can also be used in more complicated theories, e.g., in the case of 
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magnetic monopoles in non-Abelian theories ^7\. Therefore it should be possible to use these 
techniques to carry out similar studies in those systems. 

In the past, the study of dualities has often been restricted to spin models, low dimensions, 
or supersymmetric theories. The reason for this has mostly been that only in these cases 
has one been able to carry out controllable analytic calculations to test and make use of the 
dualities. Our results show how the study of dualities can in principle be extended towards 
more realistic systems, by invoking numerical techniques to corroborate analytic arguments 
based on very general principles only, such as symmetries. The same numerical techniques 
could also be used to determine "experimentally" the mappings between the parameter sets. 
It should be mentioned that many other types of numerical avenues towards duality have, of 
course, been explored and extensively tested in the case of 4-dimensional Yang-Mills theory 
(see, e.g., Ref. HD and references therein). 

To conclude, let us finally recall that apart from the theoretical issues on which we have con- 
centrated in this paper, the critical properties of three-dimensional SQED have also physical 
significance of their own, given that this theory is the effective theory of superconductiv- 
ity ^Zj and of liquid crystals [HH] in certain regimes, as well as of four-dimensional scalar 
electrodynamics at high temperatures j59 | I45 | IHTll . Particularly in the first two of these cases, 
the substantial corrections to scaling that we have observed for the photon correlation length 
(also called the penetration depth), may also have phenomenological significance jS^! given 
that it is difficult to probe very precisely the extreme vicinity of the transition point and that 
experimental samples are necessarily fairly restricted in size. 
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